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Abstract. Averaged pulse profiles of three nearby pulsars: B1929+10, J0437-4715 and B0950+08 exhibit unusual 'double 
, notches' . These W-like looking features consist of two adjacent V-shaped dips that approach each other at increasing observation 

frequency v obs roughly at a rate A oc v~jj s /2 , where A is the separation between the notches' minima. We show that basic properties 
' of the notches, namely their W-like look and the rate of their converging can be understood within a narrow class of models of 

coherent radio emission from pulsars: the free electron maser models based on coherent inverse Compton scattering of parallel 
oscillations of ambient electric field. The observed properties of the pulsars imply that the Fourier spectrum of the wiggler- 
like oscillations is narrow and that the broad-band character of the radio emission reflects the width of the electron energy 
distribution. Such a model provides a natural explanation for the frequency-independent separation between the main pulse and 
interpulse of B0950+08 as well as for the lack of radius to frequency mapping in the conal-like emission of J0437-4715. The 
frequency behaviour of the main pulse in the profile of the first radio magnetar XTE Jl 8 10- 197 can also be explained within 
J—, this model. 
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' 1. Introduction 2001). This is the opposite shift direction than expected for out- 

•rH . ward emission from a rotating magnetosphere (Blaskiewicz, 

K> .Double notches are pairs of adjacent V-shaped dips ob- ^ ■ e inm\ i \ tu. u u • t ,u 

/\ F J Cordes & Wassermann 1991). 3c) The behaviour of the po- 

!_i served in the averaged pulse profiles of three nearby pulsars: ... , . . „ ,. t . f ., 

53 ' nor. nimn m rn , 1 • c r> , nm m a in /mc s i tlon angle is not very different from the predictions of the 
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PSR B1929+10 (Rankin & Rathnasree 1997), J0437-4715 
(Navarro et al. 1997, hereafter NMSKB) and B0950+08 



rotating vector model (RVM, Radhakrishnan & Cooke 1969), 

which is puzzling because of 1) and 2), The properties 1) and 

(McLaughlin & Rankin 2004, hereafter MR04). The most strik- ~, ftl ° , . , ,. . . ■■ tf 

v b ' 2) suggest that the pedestal radio emission may originate from 

ing property of the notches is that they have the intriguing , .. . . , , , . , ... t , 

b r r 1 J b b locations in the magnetosphere that have little to do with the 

"W" look: both notches have similar (if not identical) width . , , •• • T ui j « u r 

v ' standard radio emission region. The closed held line region or 

W which is approximately equal to their separation A. The . , , .... , , . . . n 

3 n r extremely large altitudes, comparable to, or larger than, Ki c are 

separation decreases for increasing observation frequency v obs ^ exduded In the case of J04 37-4715 the notches are lo- 

(NMSKB). The notches do not affect the degree of polarization ... e • , , . u 

. . . cated in a wing of a seemingly conal component. However, the 

nor its position ang e. conal-like components in J0437-4715 also exhibit some spe- 

Remarkably, the notches reside in weak and extended cM propertieSi eg they do not follow the radius to freqU ency 

'pedestal' emission that exhibits many unusual properties: 1) mapping Infer ences of this paper refer to this special pedestal 

it is observed far from the main pulse (MP), 2) covers long and conal emission and shouM probably fee nQt extended t0 all 

intervals of pulse longitude, and 3) shows three strange polar- ]mown emission componen ts, eg. the core components, 
ization characteristics: 3a) It is very strongly linearly polarized 

(£100% for B1929+10). 3b) The S-shaped position angle curve Existing models of double notches interpret them as a dou- 
fitted to the pedestal emission only (with the MP excluded) is ble eclipse of an extended emission region by a single ab- 
shifted leftward with respect to the MP (eg. Everett & Weisberg sorber. The doubleness of such eclipse is caused by combined 
effects of differential (altitude-dependent) aberration and prop- 
Send offprint requests to: J. Dyks agation time delays within the spatially extended emission 
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region (Wright 2004). The absorber/eclipter may corotate in 
outer parts of pulsar magnetosphere (Wright) or remain sta- 
tionary at the center of the magnetosphere (Dyks et al. 2005, 
hereafter DFSRZ). 

The models based on the differential special relativistic ef- 
fects suffer from two main problems: 1) They provide no ob- 
vious reason for the W shape of the notches. Our preliminary, 
simplified calculations of pulse profiles for one version of such 
model (with the pulsar as the eclipter, DFSRZ) failed to pro- 
duce the W shape. 2) The large radial extent of the emission 
region in these models and the strong caustic effects associ- 
ated with the mechanism of the double eclipse should lead to 
rapid, complicated variations of position angle and strong de- 
polarization (Dyks et al. 2004a). None of them is observed (in 
J0437-4715 they must be caused by interaction of orthogonal 
polarization modes.) The simple property that W - A is an ex- 
tremely strong constraint on any physical model of the notches, 
and our main goal was to devise one that fulfills this require- 
ment. In Section 2 we present the observed properties of double 
notches of B1929+10 and compare it to the other pulsars. In 
Section 3 we describe a general principle of our model and the 
numerical code used to calculate modelled shapes of notches. 
Section 4 describes how the look of modelled notches depends 
on various model parameters. Sec. 5 furnishes our model with 
physics and compares it with the observations. 

2. The double notches of PSR B1 929+10 

Pulsar B1929+10 is a very useful object for studying the dou- 
ble notches because: 1) its notches are located far from the 
other strong emission components (MP, interpulse) and 2) the 
emission is strongly linearly polarized (5 100% Rankin & 
Rathnasree 1997). Both properties ensure that the pedestal 
emission within and around the notches is not contaminated by 
contributions from different emission regions. In this section 
we report new observations of B 1929+ 10 performed by one of 
us (JR) at Arecibo Observatory at frequencies 327 MHz, 1.17 
GHz, and 1.5 GHz with the respective bandwidths of 25, 100 
and 200 MHz. 

Fig. la presents the averaged pulse profile of B1929+10 
at 327 MHz. The double notches are visible at cf> = 103.5°. 
It will become clear below that a particularly constraining pa- 
rameter for theoretical models is the depth D n of the notches. 
Unfortunately, this quantity is also the most difficult to deter- 
mine from observations because it depends on the amount of 
flux received from the pulsar at the lowest (dimmest) point 
of its pulse profile. This 'unpulsed emission' was detected in 
B1929+10 via phase-resolved interferometric observations at 
408 MHz (Perry & Lyne 1985, hereafter PL85). 

To account for the depth measurement problem Fig. 1 a rep- 
resents the profile in a few different ways: The middle line 
(thin) is represented in the usual way, with the zero of the 
y-axis at the level 7 m i n of the lowest place in the pulse pro- 
file, which we take as an average of data within the phase 
range (-90, -80). For the top line we assumed that 7 m i n = 
(1.5 • 10~ 3 ± 1.5 • 10~ 4 ) 7 max , ie. we take the same ratio of 
/min/W at 327 MHz as observed by PL85 at 408 MHz. The 
error bar at <p = 90° marks the 2cr error of the baseline 
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Fig. 1. Averaged pulse profiles of B 1929+ 10 (all lines present the total 
flux), a) Three representations of the 327-MHz profile that show how 
the baseline level affects the depth £> n of the notches. The top line 
has the minimum recorded flux 7 min set at 0.0015/ max which results in 
D n =a 20%. The notches are much deeper (D„ =s 37%) if the profile is 
shown in normal way (thin middle line with / m j n = 0) . For both these 
profiles 7 max = 250. The profile at the bottom has / max = 1 and / m j n = 
0.0015/ max . b) Double notches at three different frequencies v obs . The 
horizontal bar below the notches at 1.17 GHz marks the separation 
predicted by eq. (4) for A 32 7mh z = 5.36°. The phase alignment is such 
that the MP maxima at 1.17 and 1.5 GHz lag the MP maximum at 327 
MHz by 0.8°. Different vertical scales and baseline levels were used 
at different v obs for viewing purposes. 

level from PL85. The actual error may be larger if / m i n //max 
is strongly frequency dependent. The top two lines are mul- 
tiplied by a normalization factor that sets the maximum ob- 
served flux 7 max at 250 (beyond the figure) to reveal the shape 
of the weak pedestal emission. For the third, bottom pulse pro- 
file Anax - 1 and/ m i n // max - 1.5-10 3 . The upper two curves il- 
lustrate how strongly the derived depth of the notches depends 
on the baseline level: the tiny shift of the baseline from zero 
up to 1.5 ■ 10~ 3 / max decreases D n from nearly 40% (middle 
curve) down to 20% (top). Interestingly, the flux of the pedestal 
emission increases roughly linearly with time (angle) near the 
notches. 

The 327-MHz pulse profile consists of two observations 
(12238 and 18835 single pulses long). A comparison of av- 
eraged pulse profiles for these two pulse sequences suggests 
that the relative intensity of the various weak structures in 
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discerned in published figures that present time-aligned profiles 
at different frequencies. From fig. 2 of Kuzmin et al. (1998), 
that was obtained for DM = 3.18 pc cm 4 , one can estimate 
A^mp - 1.3°. In fig. 12 of Hankins & Rankin (2006) A0 M p - 
0.3° for DM = 3.176 pc cm 4 . Interestingly, the shift of 0.8° 
puts into perfect alignment both the notches and the maxi- 
mum of the interpulse. Given that the shape of the MP changes 
strongly between 0.4 and 1 .4 GHz, it is reasonable to assume 
that the phase location of the MP maximum in B 1929+ 10 is 
frequency dependent whereas the location of notches is fixed. 
Position of double notches definitely does not depend on fre- 
quency in J0437-4715 (NMSKB) nor in B0950+08 (MR04). 
All existing observations of notches are then consistent with 
their location in pulse profiles being frequency-independent. A 
physical explanation for this will be proposed in Section 5. 

Fig. 2 presents a zoom into the same 327-MHz notches as 
shown in Fig. la (top curve). Note that their appearance (espe- 
cially the sharpness) is very sensitive to the assumed range on 
both the vertical and horizontal axis. Panel 2a presents a lin- 
ear fit to the marked data points around the notches (circles). 
The depths D n are measured between this fit and the notches' 
minima as marked with the outer vertical lines. The result is 
D\ = (20.2±1)% for the leading notch and D\ = (21.8±l)%for 
the trailing one (when expressed in percent, the depth of a given 
notch is normalized by the value of the linear fit at the phase at 
which the notch has the minimum; 2cr errors are used through- 
out this paper). In Fig. 2b the linear fit is subtracted from the 
data in two ways: the thin line represents data - fit + 1 , whereas 
the thick line shows the ratio data/fit. In either case, the lead- 
ing notch (located in the weaker emission) is shallower than the 
trailing notch in the stronger emission. The asymmetry holds at 
higher v b s . With <p — set at the MP maximum, the notches are 
located at 0i = 100.9 ± 0.07° and <j> t = 106.25 ± 0.06°, ie. they 
are separated by 5.36 ± 0.1°. The center of W lags the MP by 
103.5 ±0.1°. 

An important feature that strongly constrains the theoretical 
models is that the maximum between the double notches (at 
the center of the 'W') nearly (but not exactly) reaches the level 
interpolated with the linear fit. The data point marked in Fig. 2a 
with the central vertical line (0 = 103.5°) is 2.6% below the fit 
level, ie. it misses the fit by ~ 10% of D a . 
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Fig. 2. A zoom into the 327 MHz notches of B 1929+ 10 (expanded 
part of the top line in Fig. la). The circled points in a) have been 
used to fit the linear variations of flux around the notches. We measure 
the depth of notches D„ along the outer vertical lines; the separation 
A is the horizontal distance between them, b) The notches with the 
linear trend removed. The thin line shows data - fit + 1 ; the thick line 
presents the ratio data/fit. The leading notch is slightly less deep than 
the trailing one. c) One of model profiles (thick solid) obtained for the 
a || v case with y = 10 plotted over the data. The modelled notches are 
some ~ 10 times shallower than the observed ones. The thin solid line 
presents the same model result after a linear rescaling of the y-axis, 
which makes the modelled notches quite similar to the data. 



the pedestal may vary in time. The values of D n derived for 
these two observations are equal to 17% and 24% (assuming 
^min/^max from PL85). It is hard to tell whether this is caused 
by scintillation or intrinsic variability. 

Fig. lb shows the frequency evolution of a part of the 
profile. The notches approach each other at 1.17 GHz (mid- 
dle curve) and seem to be merged into a single feature at 1 .5 
GHz (bottom). Note that each of the profiles in Fig. lb has 
undergone different linear transformation to fit a single plot- 
ting box (no meaning should be attached to the depth of the 
notches). The profiles were phase-aligned in such a way that 
the MP maximum of the 327-MHz profile in our Fig. lb pre- 
cedes in phase the MP maximum of the 1.4-GHz profile by 
A0mp - 0.8°. Similar misalignment of the MP maxima can be 



2.1. Notches in other pulsars 

In the millisecond pulsar J0437-4715 the emission with 
notches is stronger (~ 8% of 7 max at 438 MHz) and a bit less po- 
larized (~ 70%) than in the case of B1929+10. The notches ac- 
tually seem to be carved out of the trailing wing of a prominent 
trailing component in the main pulse (see figures in NMSKB). 
At sub-GHz frequencies the top of the trailing component is 
split into two maxima separated by a single dip/notch. This bi- 
furcated top is probably bright enough (at least at v Q b s & 400 
MHz) to be studied on a single pulse basis (see Sec. 5.2). We 
are not aware of a baseline flux measurement for this object, 
so that only the upper limit of ~ 50% for the notch depth can 
be estimated. If the unpulsed flux from this pulsar is compara- 
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ble to a few percent of 7 max then the actual D a is considerably 
smaller. 

As in the case of B 1929+ 10, the double notches in 
J0437-4715 approach each other for increasing frequency v D b s , 
though they are still separable in the high quality 1 .5-GHz pulse 
profile of NMSKB. The split at the top of the trailing com- 
ponent also gets narrower at increasing frequency and looks 
nearly like a single feature at 1 .5 GHz 

On the theoretical side, the pulse profile of J0437-4715 
may be more difficult to model because its magnetosphere is 
much smaller (P = 5.76 ms, R\ c = 27.5 • 10 6 cm) than that of 
B1929+10(P = 0.2265 s,R lc = 1080- 10 6 cm) so that the radio 
emission has a larger probability of being affected by general 
relativistic effects. 

Observations of B0950+08 performed by MR04 at 430 
MHz give D n ~ 10 - 16% and PL85 report negligible amounts 
of unpulsed emission for this object at a similar frequency. 
However, the notches of B0950+08 are located close to the MP, 
in a region probably contaminated by several types of emission, 
as the low polarization degree suggests. It is therefore not ex- 
cluded that the actual D n is larger than given above. On the 
other hand, the weak emission components in B0950+08 have 
been reported to vary on a time scale of several days (MR04), 
and it is hard to tell at what stage of this variability the obser- 
vations of PL85 were done. The notches of B0950+08 seem 
fairly blurred at high frequencies (see inset in fig. 3 in MR04). 



3. The model 

The model postulates an emission region with large extent in 
the rotational azimuth and rotational colatitude 6 (A0 » A 
and A9 » A). For simplicity, we assume the emission is pro- 
duced at a fixed altitude (see Sect. 4.4 for more comments 
on emission region). The key assumption is that on a micro- 
scopic scale the region radiates a hollow cone beam along a 
local direction of magnetic field B. A single hole in the emis- 
sion region produces a cone of reduced emission, as shown in 
Fig. 3a. The notches can be observed when the observer's line 
of sight cuts through the cone. A single absorber between an 
observer and the source of radio emission will act in a similar 
way (Fig. 3b). The part of the emission region that is "half- 
hidden" from the observer has the shape of a ring. Its crossec- 
tion is marked with 'd' in Fig. 3b. 

There are at least two basic reasons for which the elemen- 
tary emission beam can have the hollow cone shape: 
1) Electrons may be accelerated along their velocity v || a 
and the radio emission observed in the pedestal can be due to 
this acceleration (Melrose 1978; Kunzl et al. 1998; Schopper et 
al. 2002; Levinson et al. 2005). In this case the opening angle 
of the elementary beam 29 R ~ 1 /y, where y ~ 10 is the Lorentz 
factor of the electrons. More exactly the received power has a 
maximum at the angle 



where /? = v/c is the electron speed in units of the speed of 
light and the approximation holds for y » 1 . In this version of 
the model the observed separation of the double notches is: 

20r 2 0.8944 1 



sin£ V5ysin<; ysin£ 



r 



(2) 



Or 



arccos 



(1 + 24/3 2 ) 1 ' 2 - 1 



4y6 



1.1^0.4472--, (1) 
5 y y 



where £ is the viewing angle measured from the rotation axis, 
the factor 1 / sin £ takes into account the 'not a great circle' ef- 
fect (eg. Dyks, Rudak & Harding 2004b), and the latter ap- 
proximation holds with accuracy better than 10% for any angle 
52° < £ < 128°. The equation is approximate in that it does 
not take into account the angle b between the hole axis and the 
line of sight. Numerical results of the next section tell us that 
this 'hole impact angle' does not affect A as long as the notches 
have the W-like shape (see Fig. 4d). Eq. (2) is valid for the case 
of the emission region with a hole (Fig. 3a). Note that accord- 
ing to eq. (2) the observed magnitude of A ~ 2° - 6° implies 
Lorentz factors of the order of 10, which is a value estimated 
on independent grounds in pulsar models based on the parallel 
acceleration (eg. the ALAE model of Melrose 1978 or CICS 
model of Schopper et al. 2002). 

2) The elementary beam has the hollow cone shape also 
when high-energy electrons (y » 10) are accelerated per- 
pendicularly (a ± v) but most of them have non-negligible 
pitch angle iff » l/y (see eg. fig. 6.5 in Rybicki & Lightman 
1979). In this case # R = i^ max , where i/f maK is the pitch angle 
for which the distribution of electrons n e (ifr) has a maximum. 
The non-zero pitch angle requires low magnetic field, such 
as is present at large altitudes, comparable to R\ c (Lyubarski 
& Petrova 1998; Malov & Machabeli 2001; Petrova 2003; 
Harding et al. 2005) as well as at lower altitudes in the millisec- 
ond pulsars (eg. J0437-4715). In spite of some positive fea- 
tures (see. Appendix A), the pitch angle version of our model 
is less successful in reproducing the data and is disfavored in 
this paper. 

3.1. The numerical code 

The numerical code to simulate the notches assumes: 1) spe- 
cific geometry and size of the emission region; 2) specific ge- 
ometry and size of the hole/fissure or absorber; 3) the structure 
of the magnetic field within the emission region, eg. dipolar, 
radial (swept out by wind) or toroidal (sweepback); 4) the ge- 
ometry of the elementary emission beam; in the case v || a it is 
calculated directly from the classical electrodynamics formu- 
lae (eg. eq. 4.101 in Rybicki & Lightman 1979) for a selected 
value of y (we do not integrate over the electron-energy distri- 
bution). In the case of v 1 a we assume the shape of the pitch 
angle distribution (eg. n e oc iff e\p[-iff 2 /i// mix 2 ], as in Epstein 
1973), and integrate the single electron 'pencil' beams of radi- 
ation over the pitch angle distribution to get the "elementary" 
emission beam for a fixed value of y. Next, the emission region 
is divided into a large number of tiny fragments, each of which 
is equipped with its own elementary beam. An observer located 
at a viewing angle £ is selected and a lightcurve is calculated by 
integration over the source for each pulse longitude (the place 
where the observer's line of sight enters the elementary beam is 
determined individually for each fragment of the source). The 
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Fig. 3. The principle of our model. 
The plane of rotational equator for a 
pulsar with dipole inclination a = 
90° chosen for simplicity is shown. 
The thick solid arc is the crossection 
through a two-dimensional emission re- 
gion. Crossections through several hol- 
low cone elementary beams of radio 
emission are indicated. The observer is 
located near the plane of the figure. The 
dotted curves mark the dipolar mag- 
netic field lines. In a the notches form 
because of the hole in the emission re- 
gion. In b they are created by a single 
absorber above the region (marked with 
bullet). The small letters help to asso- 
ciate different points of emission region 
with pulse phase at which the radiation 
emitted from them is observed. 



symmetry axis of the elementary beams is assumed to point 
along the local magnetic field B in the corotating frame. The 
aberration and propagation time delays are neglected because 
of the fixed emission altitude. 



To integrate the received flux over the electron-energy dis- 
tribution one needs to know how the coherent emissivity e at a 
specific frequency v (, s depends on the electron Lorentz factor 
y. This depends on the radio emission/coherency mechanism 
that actually works in pulsar magnetosphere. Even for a spe- 
cific emission mechanism, the dependence of e(y) at a fixed v b s 
may be completely different for different values of v bs- Let us 
take the noncoherent curvature radiation (that certainly cannot 
be responsible for the emission with notches) as an example. 
Well below the maximum in the curvature spectrum (v «: vcr) 
the emissivity does not depend on y, whereas above the max- 
imum (in the region of the exponential cut-off) it becomes ex- 
tremely sensitive to it: e C R K exp[-consf • y~ 3 ] • y~ 1/2 . For 
a weak dependence of e(y) the integration over the electron 
energy distribution could significantly blur the notches. It is 
therefore very important that in the model of the parallel accel- 
eration maser based on the coherent scattering of the 'wiggler- 
like' field, the emissivity is tightly focused around v Q b s ~ y 2 v W gi 
(see left panel of fig. 4 in Schopper et al. 2002), where v wg i is 
the frequency of oscillations of the ambient electric field. 



The angular distance ft, of the hole from the rotation axis 
and the dipole inclination a have little effect on the shape of 
notches, except from blowing them up by the factor 1/ sin ft. 
Therefore, in all cases presented in this paper we assume that 
ft = a = 90°. The emission region in all simulations of notches 
is a fragment of a sphere that extends significantly both in (p and 
direction (A0 = AO = 1O0r). 

4. Numerical results 

We have modelled numerous configurations with different hole 
geometries, absorber, B-field structure and with various pitch 
angle distributions n e (i//). In general, results are sensitive to 
bulk geometry /topology of the B-field and emission/absorption 
region. In the pitch angle case they also depend on the form of 
function « e W- 

4.1. The parallel acceleration case - a hole in the 
dipolar B -field 

Fig. 4 presents results obtained for a circular hole in the dipolar 
magnetic field (the case v || a with y = 10). In 4a the hole is 
centered at the dipole axis (6 = 0), and its angular radius ph 
increases from 0.1 to 3.3#r (top to bottom). The observer was 
located at a 'hole impact angle' b = ft - C, = 90° - £ = 0, ie. the 
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line of sight sweeps through the hole's center. For ph < 0.3#r 
the notches have the 'W'-like shape similar to the observed 
one, with the flux at the center of the 'W at nearly the same 
level as beyond the notches. An important difference between 
these cases and observations is that the depth of the modeled 
notches does not exceed a few percent, in comparison with a 
few tens of percent tentatively derived from observations. This 
is a major concern for the present version of the model and will 
be discussed below. For a larger ph the central flux drops down 
until a single wide "notch" with an initially flat bottom appears 
atp h ~0.7<9r. 

In 4b the top four curves from 4a are replotted with a 
stretched y-axis. This shows how sensitive the appearance of 
the notches is to the selected viewing method. The strong 
stretch of the y-axis produces the fall-off of flux for the increas- 
ing |0| . This occurs because the dipolar magnetic field lines are 
more spreaded at larger magnetic colatitudes 6 m measured from 
the dipole axis (so that a smaller number of lines are pointing 
towards a unit solid angle at increasing 9 m whereas the emis- 
sivity of the emission region is assumed to be uniform per unit 
surface). Note that a similar stretch of the y-axis may be un- 
knowlingly applied to plots that show the observed pulse pro- 
files if the baseline level has been overestimated (Fig. la). 

Fig. 4c shows the effect of a non-axisymmetric location of 
the hole: the dipole axis has been rotated by the angle 6 = 30° 
with respect to the line connecting the centers of the hole and 
star (see Fig. 3a). According to the basics of dipolar field ge- 
ometry, the notches moved to <p « 0.56. The displacement of 
the hole off the dipole axis results in moderate asymmetry of 
the notches' shape. The approximately linear decrease of flux 
around the notches is caused by the increasing spread of B-field 
lines. 

Fig. 4d presents the shape of notches for different viewing 
angles £ = O^-b and the other parameters (p^, y) the same as in 
the second-from-top case in panel c. For b < O.30R the 'W'-like 
notches can be observed. For larger b their shape evolves into 
a single dip. Interestingly, the separation between the notches 
practically does not depend on b as long as b < 0.70 R . The 
separation starts to decrease only for b 5 6 R , ie. when the flux 
at the W center becomes nearly as low as at the notches' min- 
ima. This property allows us to ignore the ratio of b/6 R in eq. 2 
as well as in the discussion of the frequency dependence of A 
in Section 5 below (note that b/6 R can change with v t>s for a 
single object viewed by a fixed observer). 

For the increasing y — 10, 20, 40 the modelled notches ap- 
proach each other, the flux at the W center decreases, and they 
finally merge into a single feature (solid lines in panels e and 
f of Fig. 4). This resembles the behaviour observed at differ- 
ent frequencies (Fig. lb). It can be directly interpreted within 
our model if radiation at higher frequencies is mainly gener- 
ated by electrons with larger Lorentz factor y. The depth of the 
modelled feature increases when the notches merge. This does 
not seem to happen in the data (Fig. lb) although the baseline 
level at the 1.17 and 1.5 GHz is not known. As we show in 
the next section, however, the interpretation is supported if our 
geometric idea is supplied with the physics of parallel acceler- 
ation maser. 



As we show in Appendix A one can increase the depth of 
notches up to ~ 20% in the pitch angle scenario by varying the 
function n e (i/S). However, this is usually associated with con- 
siderable deformation of their shape. Modelled notches with a 
depth and shape reasonably similar to the observed one can be 
obtained only after application of artificial "linear rescaling" of 
model lightcurves. Fig. 2c shows the result of such transfor- 
mation for the parallel acceleration case. The thick line with 
dots shows the observed 327-MHz notches of B1929+10 with 
the baseline level corrected according to the 410-MHz mea- 
surement of PL85. The thick solid line with the very shallow 
notches is the second-from-top result from Fig. 4c after a mir- 
ror reflection and horizontal shift has been applied to it (with 
no rescaling in any axis). The thin solid line in Fig. 2c is the 
same model result after a linear rescaling of y axis. This exer- 
cise is to check whether the shape of the rescaled notches in 
the a || v case bears any resemblance to the observed ones. The 
agreement is moderate: the outer sides of the modelled notches 
are less steep than the observed one. 



4.2. Geometry of the hole/absorber 

The geometry of the hole or absorber is the next important fac- 
tor that affects the shape and depth of the notches. We have 
considered a few cases that produce deeper notches. However, 
this has always been accompanied by the deformation of their 
shape. An obvious way to make the notches deeper is to ex- 
tend the hole in the direction perpendicular to the rotational 
azimuth. This, however, strongly decreases the flux at the W 
center (Fig. 5, thin line). Much the same result is obtained when 
an elongated absorber (dense plasma stream?) is placed above 
the emission region. A desperate way to avoid the decreased 
central flux is to assume an opaque thin wall/fin protruding 
from the emission region upward. Such configuration can pro- 
duce very deep notches with the central flux unaffected (if the 
fin is thin). However, the notches have a shape different than 
observed (A «c W) and become extremely asymmetric if the 
fin is located asymmetrically in the dipolar magnetic field (off 
the dipole axis). Though the fin cannot be responsible for the 
observed notches (and it would be hard to justify its origin) 
this case shows that a geometric configuration that gives deep 
notches with the unaffected central flux is possible. 



4.3. Magnetic field geometry 

In addition to the dipolar B we have considered the radial and 
toroidal magnetic field. The radial field has similar geometry 
to the dipolar field near the dipole axis and it produced similar 
results. The toroidal field (B x = cos (f>,B y = sincp, B z = 0) was 
considered as the idealized model of the swept-back, near-/?i c 
magnetic field. The direction of such B at a fixed azimuth is 
independent of 6 so that all parts of the extended source that are 
located at the meridian selected by the line of sight contribute 
equally to the observed flux. This tends to decrease D n . 
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Fig. 4. Notches modelled for a hole in 
the dipolar B-field for the case a || v. 
a) The effect of increasing hole radius 
p h for 7 = 10, b = 90° - f = and 
6 = 0, ie. with the hole centered at the 
dipole axis (top to bottom), b) Top four 
curves from a) replotted with a rescaled 
vertical axis, c) Same as b) but with 
the hole located off-axis (6 = 30°). d) 
The effect of increasing 'impact' angle 
b for p h = 0.22(9 R and y = 10. e) The 
effect of increasing y for ph = 0.57° 
and b = (solid curves). The lines 
for y = 10 and 20 were linearly trans- 
formed and replotted as the dotted and 
dashed curves, respectively, f) Same as 
e) but for b = 0.57°. Unlike in many fig- 
ures that show observed pulse profiles, 



in all cases shown with solid lines the 
zero of y-axes corresponds to 'no flux'. 



4.4. Emitter's geometry 

Our present choice of the emitter was dictated by two facts: 
1) the weak pedestal emission that contains the notches cov- 
ers a very large range of pulse phase. This implies that the 
angular extent of the emitter must be large in (/>, and conse- 
quently in 8. A negligible extent in would make the detection 
of the emission less probable, whereas the extended emission 
components seem to be quite common among the nearest and 
brightest pulsars. From the inspection of profiles in the EPN 
pulsar data base one can learn that within the distance of the 
furthest pulsar with notches (B1929+10 at 0.36 kpc) only 50% 
of pulsars in the ATNF catalogue (Manchester et al. 2005) has 
the radio emission limited to a narrow range of pulse phase 
(few tens of degrees), and some objects have really wild pulse 
profiles (eg. J2124-3358, Manchester & Han 2004). 2) The 
PA curve for the pedestal emission closely follows the curve 
of Komesaroff (1970) which probably implies that the source 
is not very extended in the radial direction (Ar <K Ri c ). The 
fixed emission altitude that we have assumed can be consid- 
ered a good approximation for any region with the radial ex- 
tent that is negligible in comparison with other relevant length 
scales. These include the scale p crv A on which the magnetic 
field changes direction by angle A, where p crv is the radius of 
curvature of magnetic field lines. The other scale is the one on 
which the electron plasma frequency drops significantly (for 
a power-law decrease of plasma density it is of the order of 
the emission height r). Both these length scales are compara- 
ble. Strong, fixed frequency radio emission from a 'too large' 
range of altitudes would probably make the notches less dis- 
tinct. The only emitter that we have managed to consider so far 
(ie. the part of spherical surface shown in Fig. 3) surely does 
not exhaust all the possible configurations that can produce the 
extended pedestal and have a reasonable detection probability. 
While simulations for other emitters are being considered, we 




0.5 r i i i - 

-20 -10 

<Pl°] 

Fig. 5. Modelled double notches for a meridionally extended fissure 
(thin line) and for a meridionally extended fin that protrudes vertically 
upward from the emission region (thick line). The notches have large 
depth but their shapes do not resemble the observed one. 



turn to the observational consequences of physical aspects of 
the model. 

5. Model versus observations 

Based on a large number of numerical results, we find that a few 
key features of model notches agree well with the observations: 
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1. Whenever the sightline cuts through the cone of reduced 
emission the calculated profiles tend to have two min- 
ima/notches; 

2. Both notches have similar depth; 

3. The width of the notches tends to be equal to their separa- 
tion. 

The list can be considerably extended if our simple idea is sup- 
plemented with the physics of a specific version of a free elec- 
tron maser: 

4. If the emission at a higher frequency v obs is mainly gener- 
ated by electrons with larger Lorentz factors, then the el- 
ementary emission beam is narrower at larger v b s which 
directly implies that the separation of notches is smaller, 
as is indeed observed. In the case of a || v, the elementary 
beam becomes narrower simply because of eq. (1). 
According to some models of pulsars (eg. Melrose 1978; 
Schopper et al. 2002; Levinson et al. 2005) the ambient 
electric field and plasma density within emission region(s) 
in pulsar magnetosphere tend to oscillate with a frequency 
y wg i that can result either from global magnetospheric elec- 
trodynamics (eg. Sturrock 1971; Levinson et al. 2005) or 
can be locally excited by streams of energetic electrons 
penetrating the ambient plasma (eg. Schopper et al. 2002). 
Hereafter we call this frequency a 'wiggler' frequency by 
a rough analogy with the laboratory free electron laser 
(FEL). 1 In the free electron maser models of coherent ra- 
dio emission from pulsars the observed radio emission can 
be considered as the Compton scattered (and blueshifted) 
wiggler frequency v wg i: 

Vobs - 7 2v wgi- (3) 

What are the spectra of y and v wg i is a difficult question, 
but one simple and seemingly natural possibility is that 
the electron energy distribution is much broader than the 
Fourier spectrum of the wiggler oscillations, ie. v wg i ~ 
const in eq. (3) for a given object. In this case different v b S 
can be associated with the inverse-Compton scattering of 
the same wiggler frequency v wg i by electrons with different 
Lorentz factors. This implies (from eqs. 2 and 3) that the 
separation of double structures should decrease with fre- 
quency according to 

A-jf JL (4, 
V5sin£\v obs / obs 

This dependence is in good agreement with the observed 
behaviour of the notches at different frequencies. Fig. 6 
presents A as a function of v bs as derived from the data on 
J0437-4715 (NMSKB, Jenet et al. 1998) and B1929+10 
(MR04; this work). All points except from the two dia- 
monds (J0437 at 1380 MHz) have 2cr errors marked. The 
locations of the diamonds have been estimated from fig. 1 
in Jenet et al. (1998). The vertical bars at these points have 
the ad-hoc length of ±10% of A. 

1 Unlike in the device, in the models mentioned above the wig- 
gler field is oscillating parallel to local B-field and electron velocity. 
Perpendicular wiggler oscillations have also been considered in the 
context of pulsars (see Fung & Kuijpers 2004). 



5. The model expressed by eq. (4) assumes that different ra- 
dio frequencies are generated by electrons that have differ- 
ent energy but occupy the same emission region. This is 
in clear contrast to the traditional view according to which 
different v obs originate from different altitudes with differ- 
ent electron plasma frequency, ie. are associated with vari- 
ations of v wg i in eq. (4). Our model predicts then that lo- 
cations of emission components in averaged pulse profiles 
should not change with v b s . This is in perfect agreement 
with the frequency independence of the separation between 
the interpulse (IP) and main pulse in B0950+08 (Hankins & 
Cordes 1981). The interpulse is connected with a bridge of 
low intensity emission with leading components of the MP. 
All the components likely have the same origin because 
they exhibit similar long term variability (cf. figs. 1 and 2 in 
MR04) as well as similar (low) intensities of single pulses 
(Nowakowski 2003). Therefore, it is probable that the en- 
tire stretch of emission that includes the IP and the 'notched 
bump' ahead of MP is generated by the inverse Compton 
parallel acceleration maser. Thus, we are surprised to real- 
ize that the issue of whether the MP -IP separation does ( or 
does not) depend on v G b s has nothing to do with the problem 
of whether we see one pole or two poles. However, since 
some parallel acceleration models have inherently built-in 
two-directional emission (eg. Levinson et al. 2005), it is 
still reasonable to interpret the structure formed by the IP, 
bridge and MP in B0950 within the two-directional emitter 
scenario (Dyks et al. 2005b; Fowler & Wright 1982; Cheng 
& Ruderman 1977). This speculation needs to be verified 
by detailed modelling. 

6. For the same reason (broad band radio emission caused by 
electron energy distribution) the conal-like components in 
J0437-4715 do not exhibit any sign of radius-to-frequency 
mapping over a very wide frequency range (McConnell et 
al. 1996; NMSKB). McConnell et al. (1996) report an in- 
trinsic low-frequency spectral turnover around ~ 100 MHz 
which can be associated within our model with a lower 
boundary of the electron energy distribution at y m [ n - 
6/ sin£. It is worth mentioning that below v Q b s - v wgiT^ in 
the relation A oc y^ 2 is not expected to hold. 

7. The oscillations of ambient E may be far from stable on 
long time scales. Evolution or disturbances of the oscil- 
lation frequency (broadening of the oscillation spectrum) 
would smear the notches. This is qualitatively consistent 
with the observed temporal variations of the shape of dou- 
ble notches (MR04; did Phillips 1990 see the notches? - 
see fig. 2 therein, and comments in Rankin & Rathnasree 
1997 or MR04). 

8 . If the emission is indeed caused by the parallel acceleration, 
there should be a chance (provided that the macrosopic ge- 
ometry of the emission region is suitable) to see the elemen- 
tary hollow cone in emission. There is indeed a bifurcated 
emission component seen in the same interval of emission 
that contains the double notches in J0437-4715 (see fig. 1 
in NMSKB). Bifurcated components observed in some pul- 
sars (eg., also for the millisecond PSR J1012+5307, see 
fig. 1 in Xilouris et al. 1998, or for the radio magnetar XTE 
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Fig. 6. Observed separation A of dou- 
ble notches as a function of fre- 
quency v obs for J0437-4715 (a) and 
B 1929+ 10 (c). The middle panel is 
for the trailing bifurcated component in 
J0437-4715. The error bars are 2tr for 
all points except from the two diamonds 
(J0437-4715 at 1380 MHz) which are 
equipped with bars of length ±10% of 
A. The straight lines mark the relation 
A oc v~' /2 predicted by eq. (4). The light 
circle for the lowest frequency point of 
J0437-4715 refers to the separation of 
notches in the averaged profile of lin- 
early polarized radiation L. The data 
are from NMSKB (circles), Jenet et 
al. (1998) (diamonds), MR04 (triangle) 
and this work (squares). 



J1810-197, see fig. 1 in Camilo et al. 2006) may result 
from the cut of our line of sight through emission from 
such a region. In Fig. 7 we present a modelled lightcurve 
for an elongated emission region that was thin in the az- 
imuth direction (A<f> = 0.3#r) but very elongated in colat- 
itude (A6 = 10#r). This configuration can be considered 
to be a zeroth order approximation of a trailing side of a 
ring centered at the dipole axis. One can see that the large 
meridional extent of the emission region does not blur the 
hollow cone shape of the elementary beam completely. The 
central panel of Fig. 6 shows that the separation between 
the maxima of the bifurcated component in J0437-4715 
also follow the relation A oc v^ 2 . 
9. The only known radio magnetar - XTE J18 10— 197 - ex- 
hibits radio emission that is variable on several long time- 
scales (days and months, Camilo et al. 2006, 2007) which 
makes it distinct from normal radio pulsars. The shape of its 
averaged pulse profile became extremely unstable in July 
2006 (Camilo et al. 2007). Several months before, however, 
the object exhibited emission with many features similar to 
those of the 'notched emission' of J0437-4715. Its radi- 
ation was highly polarized (89 + 5% at 8.4 GHz, Camilo 
et al. 2006), the main pulse was bifurcated and connected 
smoothly to an extended emission component at low v Q b s 
(see fig. 1 in Camilo et al. 2006, note that the phase align- 
ment of the profiles at different v Q b s is arbitrary). The fre- 
quency behaviour of the MP is unlike the normal radius 
to frequency mapping: instead of the evolution from a sin- 
gle component at high v b s to the well separated two conal 
components at low v b s the MP becomes broader but re- 
tains its basic shape through the very wide frequency band 
of 0.7-42 GHz. We argue here that the frequency evolution 
results from the broad electron energy spectrum (eq. 4 with 
v W gi ~ const). In Fig. 8 we show that the separation be- 
tween the maxima in the MP approach each other at a rate 
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Fig. 7. Modelled pulse profile for an emission region considerably ex- 
tended in the rotational colatitude 6 (Ad » £> R , Ai/> < 9 R ). The elemen- 
tary emission beam had a shape of the hollow cone for the parallel ac- 
celeration case with y = 10. The hollowness is revealed by the ~ 10% 
dip at the top of the pulse. The result does not depend on the viewing 
angle f . 



reasonably consistent with the v~ bs ' law at v b s ^ 1 GHz. 
The vertical error bars have magnitude of 2<x. 
At sub-GHz frequencies A ceases to increase, probably be- 
cause of a 'boundary effect' caused by a low energy limit 
y m ; n in the electron energy distribution. The sub-GHz ra- 
diation may form a low energy end of the radio spectrum 
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Fig. 8. Observed separation A of the maxima in the main pulse of the 
first known radio magnetar XTE J1810— 197 as a function of observa- 
tion frequency v obs . The errors of A have magnitude of 2cr. The hori- 
zontal bars (too short to be visible for most points) show a bandwidth. 
The straight lines mark the relation A oc v~ b ' s ' 2 given by eq. (4). The 
data are from Camilo et al. (2006). 

with v b s ~ v wgi7min- The three flux-density measurements 
between 0.69 and 2.9 GHz done on MJD 53850.9 (Table 1 
of Camilo et al. 2006) are indeed consistent with a harder 
spectrum (S v oc yJJ ) than measured at higher frequencies 
(S v oc v ob ° s 5 betwee°n S 1.4 and 19 GHz) on MJD 53857. 
10. Contrary to other models of double notches (Wright 2004; 
DFSRZ) the emission region considered here does not have 
to be radially extended (it may extend moderately or not at 
all). The large depolarization and extremely complicated 
position-angle curve are characteristic features of regions 
with large radial extension (Dyks et al. 2004a) and are not 
observed in the pedestal emission. The observed polariza- 
tion properties of the pedestal emission thus seem to be 
more consistent with the present model. 



Had it been possible to observe both the double notches and the 
periodic microstructure with the timescale r wg i, one could use 
the above equation to constrain 

The width of the wiggler spectrum allowed by the model 
can be assessed as follows: If we assume that the notches get 
blurred when 1 jy changes "significantly" (say by a factor of 2), 
then the spectrum of v wg i cannot extend by a factor larger than 
4 (since 1/y oc v^ 2 ). This rough estimate applies only to the 
part of magnetosphere that produces detectable radio emission. 

For B1929+10 we observe A 5.3° at 327 MHz, ie. y = 
9.7/ sin £ and r wg i = 0.28 /is sin _2 f which requires sin- 
gle pulse observations with sub-microsecond time resolution. 
Given the weakness of the pedestal emission in B1929+10 
this observation can be unfeasible even for SKA. In the other 
pulsars, however, the emission with notches (or with double 
emission features) is much brighter (a few - 30 percent of 
Anax) and should be easier to observe. For B0950+08 one gets 
A 5.3° at 430 MHz (MR04) so that y 9.7/ sin^ and 
r wg i = 0.22 yus sin -2 f . The maxima of the bifurcated trailing 
component in the pulse profile of J0437-4715 are separated by 
some 3.9° at 438 MHz (NMSKB). This gives y 13/ sin£ and 
i" W gi -0.4/is sin" 2 £ = 1. 6 fis (sin 30°/ sin £)~ 2 . Attempts to re- 
veal the periodicity of microstructure in J0437-4715 have been 
done at two widely separated frequencies (327 MHz, Abies et 
al. 1997, and 1380 MHz, Jenet et al. 1998) with apparently 
conflicting results. The periodicity of 22 /is reported by Abies 
et al. would be consistent with eq. (5) for a very small view- 
ing angle f ~ 8° that we consider somewhat extreme. Jenet 
et al. (1998) report no microstructure periodicity down to 80 
ns but their single pulse signal could have been dominated 
by the bright core component with little contribution from the 
'notched emission' that is weak at 1380 MHz. 

The "absorber version" of our model allows for a range of 
possible microstructure periodicity, because the counterpart of 
eq. (5) becomes additionally dependent on geometric param- 
eters of the system (eg. the distance of the absorber from the 
emission region, the curvature of the region, etc). Whereas the 
v^ 2 dependence can appear in this case under some conditions 
(eg. in radial B-field), a specific A translates in this scenario 
into a larger opening angle Or and smaller y (see Fig. 3). The 
expected r wg i are therefore smaller than in the case shown in 
Fig. 3a. 



5.1. Predicted periodicity of microstructure 

The microstructure observed in pulsars exhibits typical 
timescales or quasi-periodicities (eg., Cordes et al. 1990; Lange 
et al. 1998) that used to be interpreted within the FEL type 
models of pulsars as a direct result of the wiggler oscillations 
and are used to estimate the value of the Lorentz factor in 
these models. Since in our model y can be independently es- 
timated from the notches' separation, the model can be verified 
by checking whether the typical timescale of microstructure is 
equal to 

i r 2 Q-8 ... 

1"wgl = - - — . 2 • {•>) 

Vwgi v b s Vobs A z sin 4 



5.2. Single pulse visibility of double notches and 
emission cones 

Our numerical simulations assumed that the radio emissivity 
was uniform and steady throughout the entire emission region, 
because we were modelling double features in averaged pulse 
profiles. A natural question is whether the double notches can 
be observed in single pulse emission. A related question is 
whether the hollow cone shape of the elementary radio beam 
can be directly observed as pairs of emission features in the in- 
stantanenous (single pulse) radio emission. The simple answer 
to these questions is 'no' - the notches that are pronounced 
in an average profile should not be seen as an absorption fea- 
ture in the single pulse data. The hollowness of the emission 
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cone that can be evident in the averaged pulse profile (Fig. 7) 
cannot be recognized in single pulses either. The reason is that 
the instantaneous emissivity within the emission region is very 
non-uniform and variable as suggested by features observed 
on a variety of timescales shorter than P (eg. Weltevrede et 
al. 2006; Edwards et al. 2003; Johnston & Romani 2002; Cairns 
et al. 2004) and as is normally assumed in the models of drift- 
ing subpulses (Ruderman & Sutherland 1975; Wright 2003). 

The detailed discussion is deferred to Appendix B. Here we 
only mention the timescale relevant to the problem to isolate 
some limiting cases. To simplify the analysis we assume the 
equatorial geometry illustrated in Fig. 1 (a — £ — 90°). Let us 
consider a single bunch of electrons (a spark) associated with 
a single hollow cone of radio emission of opening angle 20 R . 
Let us initially assume that the bunch is 'frozen', ie. it does not 
move in the corotating frame and its emission cone does not 
evolve. If the leading side of its cone is directed towards the 
observer at some moment t\ then a considerable period of time 



Af m , - 



20 R PA 



Q 



2tt 



7 • 10" 



P A 

5ms 5° 



(6) 



must pass before the rotation of the magnetosphere directs the 
trailing side of the cone towards the observer, at moment t t . To 
detect both sides of the cone the fictious frozen bunch would 
have to survive for At Z Af rot ~ 7 10~ 5 - 3 10~ 3 s, where 
the range corresponds to P 5 - 250 ms observed among the 
pulsars with notches. The corresponding light travel distance is 
Ar rot = cAf rot ~ (2 — 105) • 10 6 cm, ie. of the order of, or much 
larger than the neutron star radius R N $. In Appendix B we dis- 
cuss two types of limitations that make the observability of the 
notches in single pulses improbable/impossible: one is purely 
geometrical and associated with the non-uniform illumination 
of the radially thin emission region (slowly drifting bright spots 
that do not evolve on timescales shorter than Af rot nor do they 
move relativistically in any direction). The other one is due to 
special relativistic kinematics and refers to a realistic situation 
of fast moving (outflowing) bunches of electrons that evolve on 
a timescale At <sc Af rot that is too short for a bunch to expose 
both sides of its elementary emission cone to the observer. 



6. Conclusions/Discussion 

We conclude that a remarkable number of peculiar observa- 
tional effects can be understood within a radio coherency model 
that is based on inverse-Compton scattering of a single 'wig- 
gler' frequency (or a narrow band of wiggler frequencies) by a 
broad energy distribution of electrons. The model explains the 
'W'-like shape of the notches, the bifurcated emission compo- 
nents, the convergence of these features at the rate v~ b ' s /2 , the 
frequency independence of the separation between the IP and 
MP in B0950+08 as well as the lack of radius-to-frequency 
mapping in J0437-4715. The model performs reasonably well 
for the only known radio magnetar XTE J 1 8 1 0- 1 97, which im- 
plies that the same mechanism of coherent radio emission is 
operating in objects with so wildly different surface magnetic 
fields as the magnetars and millisecond pulsars. This finding 
is consistent with the linear acceleration origin of radio emis- 



sion, because the strength of B-field is largely irrelevant for this 
mechanism (Rowe 1995). 

There are many unsolved puzzles that remain and that need 
to be addressed in future. These are the depth of the notches, 
the macroscopic geometry of the emitter and hole/absorber, the 
nature of the hole, etc. They are closely interrelated so it may 
be worthwhile to address all of them simultaneously rather than 
treat them separately. The large depth of double notches could 
easily be produced if the coherent radio emission had occured 
at two small angles (~l/y) with respect to the plane of B-field 
lines. When smeared along the B-field lines this would result 
in a kind of two-planar emission geometry. The maser based 
on curvature drift (Zheleznyakov & Shaposhnikov 1979; Luo 
& Melrose 1992, hereafter LM92) offers this kind of emis- 
sion geometry (see the bottom panel of fig. lb in LM92) but 
it predicts much weaker frequency dependence of A than ob- 
served (A oc v^ 3 according to eq. (9) of LM92). In the case of 
models based on parallel acceleration, part of the hollow cone 
beam near the B-field line plane would have to get absorbed to 
obtain the two-planar emission geometry. An interesting ques- 
tion is whether the emission is outward or perhaps inward, as 
suggested by the leftward shift of the position-angle curve. It 
is worth emphasizing that two-directional emission is inherent 
in some of the parallel acceleration models (eg. Levinson et 
al. 2005), although bulk of a magnetosphere of fast rotating 
pulsars (J0437-4715) has been estimated opaque for the back- 
ward emission (Luo & Melrose 2006). An important issue is 
the altitude of the emission. The typical radio frequency ex- 
pected in free electron maser models is proportional to electron 
plasma frequency and is too high in comparison with the ob- 
served ~ 100 MHz emission (Kunzl et al. 1998; Melrose 2000). 
This favors large emission distance (r ~ R\ c ) which is con- 
sistent with the far-from-MP location of the pedestal emission 
(provided that B0950+08 or B 1929+ 10 have large dipole in- 
clinations, which seems probable). The strength of the guiding 
magnetic field required to stabilize parallel Langmuir waves 
is too low (10 3 G in Schopper et al. 2002) to constrain pos- 
sible emission site. On the other hand, the nearly RVM shape 
of the position-angle curve (Rankin & Rathnasree 1997) sug- 
gests low emission altitudes and these would require radiation 
from the closed field line region. The region is a place of copi- 
ous pair production according to the outer gap model (Cheng 
et al. 1986; Hirotani et al. 2003; Takata et al. 2006; Wang et 
al. 1998, see fig. 3 therein), which is quite successful in re- 
producing gamma-ray pulse profiles of pulsars (eg. Romani & 
Yadigaroglu 1995; Dyks & Rudak 2003). The details of the 
macroscopic geometry of the system remain a puzzle, and we 
emphasize that the configurations shown in Fig. 3 may be very 
far from reality. 

A big unsolved issue is what is the relation of the radio 
emission considered in this paper to the more normal (?) emis- 
sion that can be classified within the scenario of core and conal 
beams that exhibit nulling, drifting as well as the radius to 
frequency mapping. Can the RFM-exhibiting conals be inter- 
preted within the same 'parallel FEL model, but dominated by 
variations of v wg i in eq. (3)? What is the origin of core emis- 
sion? 
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220 230 240 250 

pulse phase [deg] 

Fig. 9. A fragment of pulse profile of B0950+08 that shows the be- 
haviour of double notches at increasing frequency (top to bottom). 
The data are from the European Pulsar Network base (Gould & Lyne 
1998). Pairs of vertical bars above each profile present the relation 
A = 6.5°(408 MHz/v obs )" 2 (eq. 4). 

The model proposed here can be tested observationally 
by searching for microstructure periodicities in the emission 
with notches or with bifurcated components (eq. 5). The ex- 
pected timescale of the microstructure is of the order of 1 
/is. Multifrequency observations of other objects at a high 
signal-to-noise ratio can provide further support for the rela- 
tion A oc v'V 1 . Since there is some evidence of temporal evo- 
lution of double notches, their separation at a fixed frequency 
may in principle vary in time. Therefore, simultaneous mul- 
tifrequency observations would be most preferred, although 
even the non-simultaneous data in Figs. 6 and 8 proved suc- 
cessful in revealing the A(v Q b s ) relation. The millisecond pul- 
sar J1012+5307 has a bifurcated component and is a very good 
candidate (see. fig. 5 in Kramer et al. 1999). The double notches 
of B0950+08 present a difficult observing target (Fig. 9) and it 
is important to determine their multifrequency behaviour with 
better definition. The broad band approach can also tell us 
whether breaks in the relation A(v Q b s ) are associated with spec- 
tral breaks. A very important but difficult task is to precisely 
determine baseline levels for the pulsars with notches at dif- 
ferent frequencies. This would provide the depths of notches 
that strongly constrain possible geometric configurations of the 
magnetospheric emitter. 

Acknowledgements. This paper was inspired by the recent works 
of Geoffrey A. E. Wright (MNRAS 2003, 2004). JD thanks 



J. Borkowski, F. Camilo, Q. Luo, and D. Melrose for helpful 
comments. I am also grateful to R. Manchester and S. Ransom 
for the amazing pulsar data. We acknowledge the support of the 
KBN grant N203 017 31/2872 and NSF grants (JMR). Arecibo 
Observatory is operated by Cornell University under contract to 
the US National Science Foundation. Part of this research has 
made use of the data base of published pulse profiles maintained 
by the European Pulsar Network, available at: http://www.mpifr- 
bonn.mpg.de/pulsar/data/. We have also used the ATNF pulsar 
Catalogue (http://www.atnf.csiro.au/research/pulsar/psrcat). 

References 

Abies, J.G., McConnell, D., Deshpande, A.A., & Vivekanand, M. 

1997, ApJ, 475, L33 
Blaskiewicz M., Cordes J.M., Wasserman I., 1991, ApJ, 370, 643 
Cairns, I.H., Johnston, S., & Das, P. 2004, MNRAS, 353, 270 
Camilo, E, Ransom, S. M., Halpern, J. P., Reynolds, J., Helfand, D. 

J., et al., 2006, Nature, 442, 892 
Camilo, E, Cognard, I., Ransom, S. M., Halpern, J. P., Reynolds, J., et 

al., 2007, ApJ, submitted (astro-ph/06 10685) 
Cheng, A. E, & Ruderman, M. 1977, ApJ, 216, 865 
Cheng, K. S., Ho, C, & Ruderman, M. A. 1986, ApJ, 300, 500 
Cordes, J. M., Weisberg, J.M., & Hankins, T.H. 1990, AJ 100.2, 1882 
Daugherty J.K., Harding A.K., 1983, ApJ 273, 761 
Dyks, J., Frackowiak, M., Slowikowska, A., et al., 2005, ApJ 633, 

1101 (DFSRZ) 
Dyks, J., Zhang, B., & Gil, J. 2005b, ApJ, 626, L45 
Dyks J., Harding A.K., & Rudak B., 2004a, ApJ 606, 1 125 
Dyks J., Rudak B., & Harding, A. K. 2004b, ApJ 607, 939 
Dyks, J., & Rudak, B. 2003, ApJ, 598, 1201 

Edwards, R.T., Stappers, B.W., & van Leeuwen, A.G.J. 2003, A&A, 
402, 321 

Epstein, R. I., 1973, ApJ, 183, 593 

Everett J.E., Weisberg J.M., 2001, ApJ 553, 341 

Fowler, LA., & Wright, G.A.E. 1982, A&A, 109, 279 

Gould, D.M., & Lyne, A.G. 1998, MNRAS, 301, 235 

Hankins, T.H., & Cordes, J.M. 1981, ApJ, 249, 241 

Hankins, T.H., & Rankin, J.M., 2006, preprint 

Harding, A.K., Usov V., Muslimov A.G., 2005, ApJ 622, 531 

Hirotani, K., Harding, A.K., & Shibata, S. 2003, ApJ, 591, 334 

Jenet, FA., Anderson, S.B., Kaspi, V.M., et al. 1998, ApJ, 498, 365 

Johnston, S., & Romani, R.W. 2002, MNRAS, 332, 109 

Komesaroff, M.M. 1970, Nature, 225, 612 

Kramer, M., Lange, Ch., Lorimer, D.R., et al. 1999, ApJ, 526, 957 
Kunzl T., Lesch H., Jessner A., von Hoensbroech, 1998, ApJ 505, 
L139 

Kuzmin, A.D., Izvekova, V.A., Shitov, Yu.P, et al. 1998, A&AS, 127, 
355 

Lange, Ch., Kramer, M., Wielebinski, R., & Jessner, A. 1998 A&A, 
332, 111 

Levinson, A., Melrose, D., Judge, A., & Luo, Q. 2005, ApJ, 631, 456 
Luo, Q., & Melrose, D. B. 1992, MNRAS, 258, 616 
Luo, Q„ & Melrose, D. B. 2006, MNRAS, 371, 1395 
Lyubarski, Y.E., & Petrova, SA. 1998, A&A, 337, 433 
Malov, I. E, & Machabeli, G. Z. 2001, ApJ, 554, 587 
Manchester, R.N., Han, J.L. 2004, ApJ, 609, 354 
Manchester, R.N., Hobbs, G.B., Teoh, A., & Hobbs, M. 2005, Astron. 
J., 129, 1993 

McConnell, D., Abies, J.G., Bailes, M., et al. 1996, MNRAS, 280, 331 
McLaughlin, M. A., & Rankin, J. M. 2004, MNRAS, 351, 808 
(MR04) 

Melrose, D. B. 1978, ApJ, 225, 557 



J. Dyks et al.: A model for double notches and bifurcated components in radio profiles of pulsars and magnetars 



13 



Melrose, D. B. 2000, ASP Conference Series, 202, 721 

Navarro, J., Manchester, R. N, Sandhu, J. S., et al., 1997, ApJ, 486, 

1019 (NMSKB) 
Nowakowski, L. 2003, Arecibo Newsletter, 36, 8 
Perry, T.E., & Lyne, A. G. 1985, MNRAS, 212, 489 (PL85) 
Petrova, S. A. 2002, MNRAS, 336, 774 
Phillips, J. A. 1990, ApJ, 361, L57 

Radhakrishnan, V., & Cooke, D.J. 1969, Astrophys. Lett., 3, 225 

Rankin, J. M., & Rathnasree, N. 1997, J. Astrophys. Astron., 18, 91 

Romani, R.W., & Yadigaroglu, I.-A. 1995, ApJ, 438, 314 

Rowe, E.T. 1995, A&A, 296, 275 

Ruderman M.A., Sutherland PC 1975, ApJ 196, 51 

Rybicki G.P, Lightman A.P, 1979, "Radiative processes in 

Astrophysics", Wiley-Interscience, New York 
Schopper R., Ruhl H., Kunzl T.A., Lesch H., 2002, MPE Report 278, 

193 

Sturrock, PA. 1971, ApJ, 164, 529 

Takata, J., Shibata, S., Hirotani, K., & Chang, H.-K. 2006, MNRAS, 
366, 1310 

Wang, F.Y.-H., Ruderman, M., Halpern, J.P, & Zhu, T. 1998, ApJ, 
438, 314 

Weltevrede, P., Edwards, R.T., & Stappers, B.W. 2006, A&A, 445, 243 
Wright, G. A. E., 2003, MNRAS, 344, 1041 
Wright, G. A. E., 2004, MNRAS, 351, 813 

Xilouris, K. M., Kramer, M., Jessner, A., et al. 1998, ApJ, 501, 286 
Zheleznyakov, V. V., & Shaposhnikov, V.E. 1979, Aust. J. Phys., 32, 
49 

Appendix A: The pitch angle case 

The pitch angle case is discussed here for the sake of complete- 
ness, as well as in view of some positive features it has (see 
below). We begin with a short description of numerical results 
and then discuss how the model performs in confrontation with 
reality. 

As the simplest choice of the electron's pitch-angle distri- 
bution n e (iff) we take the triangular shape: 

n e (t/0 = for < ifr < i/w - Aif/ 

"eW = c in [if/ - (lf/ m . dx - Ai/OJAi/T 1 

for i/w - Ai// < \jj <\p 

max / * 1 \ 

«eW = C out [i// - (</W + At/OK-At/T 1 ) 

for t/w < tfr < i/w + Ai/' 

n e (t/0 = for i/w + Aif/< (// <n. 

where c; n and c out can take the values of or 1 and Aij/ is the 
bottom width of the distribution. For c; n = c out = 1 the value 
of n e (if/) increases linearly from at i^ max - Aif/ up to 1 at ^ max , 
then drops linearly down to zero at if/ m . dx + Aij/. For c; n = (and 
Cout = 1) only the outer part of the distribution is present with 
a sharp inner boundary at i^ max . For c out = (and c; n = 1) the 
inner part of the distribution preserves. 

Fig. A.l presents the look of notches for various types of 
the pitch-angle distribution n e W- In A. la we have assumed a 
'two-sided' distribution (c; n = c out = 1) with i^ max = 5° and 
Ai/' = O^t/'max which is an example of a general form of n e W 
that is closely confined to ^ max (Aif/ <§: ^ max ). This case differs 
from the parallel acceleration case mainly in that n e (if/) covers 
much smaller solid angle (in comparison to 9^) than the beam 
of non-negligible emission does in the v || a case. This has two 
consequences visible in Fig. A. la: 1) The notches retain the 



high flux at the W center for the hole size larger than in the 
parallel acceleration case (p n ~ 0.5#r) because no radiation is 
emitted near the symmetry axis of the pitch-angle distribution 
(iff < 4°). 2) The W-shaped notches (with the unaffected cen- 
tral flux) are much deeper and reach the ~ 20% depth observed 
for B1929+10. However, their shape becomes very distorted in 
comparison to the observed one. For p n < 0.4# R the notches 
are narrow and separated by a flat plateau so that A » W. For 
larger p n the outer sides of the notches become rounded and 
much less steep than the inner ones that form the central max- 
imum. All curves in panel a) correspond to the central viewing 
ofthehole(£> = 0)andy = 10 3 . Fig. A. lb presents resultforthe 
same pitch-angle distribution as in a) but p n is fixed to 0.55#r 
and b is variable. 

In Fig. A.lc i^ max = Ai// = 5°, ie. « e (<A0 has a width com- 
parable to tf/ mdx . As before, c m = c out = 1, ie. the distribution 
is symmetrical with respect to t// max . This shape is qualitatively 
similar to the shape of the radiation beam for the parallel ac- 
celeration case. Therefore the resulting notches look alike. The 
result is for b — 0. 

Fig. A. Id presents results for the inner pitch-angle distribu- 
tion (c out = in eq. A.l) in the case of b = and variable p n . 
Here the separation between the notches is in general different 
from their width: A > W for p h ^ 0.29 R whereas A < W for 
Ph ~ O.30R. The notches have peculiar shape. 

Fig. A.le is for the outer n e W distribution (eq. A.l with 
c; n = 0). The shape of notches is again unlike the observed one. 
The W center is flat for p n $ 0.5#r or has a dip for p n - O.750R 
and the notches look strange. 

Fig. A. If is for the pitch-angle distribution considered by 

Erber (1973): n e K i^exp [-tA^/iAinax]- ^ ts § enera l properties 
are similar to the case shown in Fig. A.lc (At/' ~ ^ max ) and 
therefore the resulting notches have similar shape, close to the 
observed one. 

The general property that is common for all the cases is 
that for very small p n the hole acts as a delta function and the 
resulting notches have the same shape as n e W- Therefore, the 
notches have sharp minima (and are very shallow) whenever 
Ph ^ i^max (topmost curve(s) in panels a, c, d, and e). For larger 
Ph the minima become oblate in all cases with the exception of 
the inner n e (i//) case (panel d). It is worth emphasizing that the 
observed notches appear much sharper than they really are on 
most plots that show full rotation period (cf. top line in Fig. la 
with exactly the same profile shown in Fig. 2a). 

A few points on confrontation of the model with the data: 

The conclusions 1 to 3 in section 5 are valid for both the 
parallel acceleration case as well as for a (quite numerous) class 
of the pitch-angle distributions that fulfill i^ max ~ At}/. 

The initial pitch angle of electron-positron pairs created 
through one-photon absorption is also inversely proportional 
to y. This is because photons propagating at small angles 
if/ relative to B need to have larger energy to produce pairs 
[(e/mc 2 )(B/BQ)smif/ * 2/15, Ruderman & Sutherland 1975] 
and the pair components tend to share the energy e of a par- 
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Fig. A.l. Modelled double notches cal- 
culated for various pitch-angle distribu- 
tions with a fixed y = 10 3 and i^ max = 
5°. A hole in the emission region is as- 
sumed. For details see text. 



ent photon equally (jmc 2 « e/2, Daugherty & Harding 1983), 
which leads to 



Afterwards ifr undergoes strong evolution but its value remains 
anticorrelated with y (see fig. 2 in Harding et al. 2005). 

The distribution of pitch angles may be considered fragile 
and susceptible to disturbances, which is consistent with the 
probable variability of the notches and pedestal emission on a 
timescale of days. 

The low \B\ required for the large i^ max is consistent with 
the observed location of the pedestal emission (far from MP 
in B1929+10 and B0950+08) as well as with the presence of 
notches in J0437-4715. The magnetar does not seem to fit into 
the picture, although we do not know how large the radio emis- 
sion altitude in XTE Jl 8 10- 197 might be. 

Since we do not see how the relation A oc y^ 2 could arise 
in the scenario of the favored pitch angle, we consider it less 
natural than the a \\ v case. 

Appendix B: Single pulse visibility of double 
features 

B.1. The quasi-steady but nonuniform source at fixed 
altitude (drifting spots) 

As can be seen in Fig. 3b the notches can be observed in full 
form in single pulse data only when the absorber is nearly- 
simultaneously (or: for a sufficiently long period of time > 
Af rot ) illuminated from various directions by a sufficiently large 
part of the emission region (extending at least between the two 
points marked with the small letters 'd' in Fig. 3b). In the 
hole case a nearly-simultaneous emission around the hole is 
required (or at least on opposite sides of the hole.) For a ran- 
dom distribution of sparks (emission spots) of angular radius 



Psprk ^ #r one can observe only a part of the notch feature 
(a shoulder or dip) or no absorption feature at all. The latter 
should frequently occur because there are a great deal of points 
on the emission region (eg. marked with 'b' in Fig. 3) that con- 
tribute unobscured radiation at the pulse phase of the notches. 
The phase at which the radiation is detected is marked also with 
'b' just above the schematic pulse profiles at the top of Fig. 3. 
In this scenario it would be natural to see no single pulse with 
the double notches. There would also be lots of single emission 
features because of contribution of points located outside of the 
plane of Fig. 3 at angular distance ~ # R from the line of sight 
traverse. One would have also seen, however, a great number 
of double subpulses corresponding to the sightline cut through 
the elementary cones, which is not observed. The model just 
discussed is thus excluded by the observed single pulse proper- 
ties of pulsar radio emission. In principle, however, such kind 
of emitter could be a part of magnetosphere that stays at a fixed 
altitude and is steadily refreshed by electrons that keep flowing 
through it. 

B.2. Bunches of electrons outflowing at relativistic 
speeds 

A much more realistic model for the basic emission units in 
pulsar magnetosphere is based on short living bunches of elec- 
trons that outflow at relativistic speeds (y ~ 10). Such sources 
are localized both in the radial and horizontal direction, but 
their radial coordinate increases at nearly the speed of light 
dr/dt = 0.995c for y = 10. The chance to observe both sides 
of the elementary emission cone from a single bunch fully de- 
pends on whether the bunch can survive for sufficiently long 
period of time, namely the time interval f t - h needed for the 
cone to expose its other side to the observer. 

The source (bunch) is now approaching the observer (and 
trying to catch up with the photons) with y ~ 10 so that its ra- 
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dial distance changes according to r ^ ct + ro. For definiteness, 
hereafter we assume that the first side of the cone is spotted 
by the observer at time to = t\ = and at the radial distance 
r a = K f o) ~ ^ns that will be neglected because r <K Ar rot for 
most periods we are interested in. Let the azimuth of the emis- 
sion direction on the observed side of the cone is (f>\ — (in the 
observer's frame OF), which is also the fixed azimuth of the 
line of sight. For a bunch that flows along the dipole axis, the 
azimuth of the other side of the cone will change according to 



t - 29 R - Q.t ■ 



flic 



A-2 



ct 



(B.l) 



where 26 R is the initial value, Qt takes into account the ro- 
tation of the dipole and r/R\ c takes into account the forward 
projection of radiation caused by the aberration effect in the 
limit r <sc r (eg. Dyks, Rudak & Harding 2004b). The ob- 
server can see the trailing side of the cone (</> t = 0) after 
time t = PA/ (An) which is two times smaller (3 • 10~ 5 - 10" 3 
s) than Af rot given by eq. (6) and corresponds to light travel 
lengthscales of (1 -50)- 10 6 cm. The scales are much longer 
than predicted/considered in some models of radio coherency, 
eg. Ar ~ a few x 10 5 cm in Schopper et al. (2002) or in the 
perpendicular acceleration model of Fung & Kuijpers (2004). 

Had, however, the bunch survived for the time f t - 1\ without 
any significant evolution (to maintain its original hollow cone 
emission beam) then the trailing side of the cone would not be 
observed ~ 5° after the leading side, but much sooner. Because 
the bunch propagates towards the observer with the speed v 
close to c, the radiation in the trailing side lags the radio waves 
from the leading side by Ar cAt - vAt cAt/(2y 2 ). The 
trailing radio waves are then detected just after the leading side 
radio waves: 



Afobs - 



2y 2 



5 ■ 10"* s 



10" 



\10/ 



(B.2) 



so that nAf b s «c A. To detect both sides of the hollow cone 
from a single bunch or particle, one would therefore need to 
look for double emission features on the timescales of nanosec- 
onds to milliseconds. The bunch, however, and the physical 
conditions around it could not evolve significantly during the 
time ? t - 1\ ~ 10~ 5 - 10- 3 s which may be impossible to satisfy. 

We therefore conclude that the direct observation of the 
elementary hollow cone in the single pulse data is extremely 
improbable, if not impossible. In any case, in the outflowing 
bunch scenario the ~ 5° opening angle of the cone would corre- 
spond (in single pulse data) to double emission spikes separated 
by ~ 0.01°. The 5° separation is visible in the averaged pulse 
profiles because the trailing notch (or the trailing maxium in 
Fig. 7) is created by radio waves emitted from roughly the same 
altitude as those associated with the leading notch/maximum, 
but emitted later by a different bunch of electrons. Therefore, 
the only sign of the double notches in the single pulse data 
will be the less frequent appearance of single-looking, narrow 
emission spikes that are normally observed in the high time 
resolution data. 



B.3. The role of B -field line curvature 

The above estimates have been done for the special case of the 
dipole axis, but the curvature of magnetic field lines can drag 
the other side of the hollow cone into the observer's view. On 
the leading side of the dipole axis, the eq. (B.l) becomes: 



h * 20r - at - ^ — 

«lc Pcrv 



(B.3) 



where p crv in the last term is the radius of curvature of mag- 
netic field lines. The second side of the hollow cone is directed 
towards the observer after time 



ct 



Perl + ■ 



(B.4) 



(4/3)(rRy.) 1 ' 2 = 9.2 



At the last open field lines p crv 
10 7 (Pr6) 1/ ' 2 where r^ = r/(10 6 cm). For the considered range 
of periods (5 - 250 ms) and r ~ 10 6 cm the resulting length- 
scales are (0.5 - 5) • 10 6 cm and become larger for increasing 
r, eg. (0.7 - 12) • 10 6 cm for r(, = 10. The lower limit of the 
scales (that refers to the 5 ms period) is still of the order of R^. 
The scale has decreased considerably only for the normal, long 
period pulsars. 

On the trailing side of the dipole axis the sign of the last 
term in eq. (B.3) is positive and the equation has no positive 
solutions for t if p a l v > 2RZ 1 , eg. near the edge of the polar cap. 
This simply means that the combined effects of rotation and 
aberration (2nd and 3rd terms in eq. B.3) are too weak to com- 
pensate for the backward bending of the B-field lines, ie. the 
field line's curvature drags the beam away from the observer's 
line of sight. A full cut through the beam is possible only when 
the observer sees the trailing side of the cone first. Appropriate 
change of signs in (B.3) gives then 



ct 



Perl - 2R- 1 ' 



(B.5) 



which gives timescales only slightly larger than on the lead- 
ing side because p crv «: R\ c /2 in the region where we make 
estimates (low altitudes, close to the last open field lines). An 
interesting solution of (B.5) corresponds to p c „ = R\ c /2 (field 
lines slightly on the trailing side of the dipole axis) for which 
t — 00. This is the case in which the beam does not rotates in 
the obsever's frame, ie. the field line curvature fully compen- 
sates aberration and rotation. The observer can see the leading 
side of the beam all the time as the bunch propagates upwards 
(caustic pile up of radiation near a fixed phase in pulse profile), 
actually up to the altitude above which our simple approxima- 
tions (B.1)-(B.5) break down. 

In the curved magnetic field lines the hollow cone emission 
caused by the parallel acceleration becomes dominated by the 
curvature emission at frequencies v cj 3 /(2np crv ) (Melrose 
1978). For the range of P = 5 - 250 ms and p crv at the polar 
cap rim the curvature radiation declines above 0.1 - 0.7 MHz 
which is well below the frequencies at which the notches are 
observed. 



